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For  a  pair  of  reaction  diffusion  equations  with  one  diffusion  coefficient 
very  large,  there  is  associated  a  reaction  diffusion  equation  coupled  with  an 
ordinary  differential  equation  (the  shadow  system)  with  nonlocal  effects  which 
has  the  property  that  it  contains  all  of  the  essential  dynamics  of  the  original 

•>. .  .  '"/■  ••  •  i  ’  i ..  .  , 


equations. 


V  ■  '  ' 

r  r  ‘ 
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I.  introduction 


Many  models  of  chemical,  biological  and  ecological  problems  involve 
systems  of  reaction-diffusions  in  a  bounded  domain  ft  with  Neumann 
bound  >  '-onditions.  Of  major  concern  is  an  understanding  of  the  mechanism 
for  the  creation  of  stable  patterns;  that  is,  stable  solutions  which  are  spatially 
dependent  (see,  for  example,  Turing  [1952],  Nicolis  and  Prigoginc  [1973]).  For 
the  understanding  of  how  stable  patterns  occur,  it  is  obviously  of  interest  to 
characterize  those  situations  for  which  stable  patterns  do  not  exist  and,  even 
more  particularly,  those  systems  for  which  the  flow  is  essentially  determined  by 
an  ordinary  differential  equation.  This  situation  was  studied  rather  extensively 
by  Conway,  Hoff  and  Smoller  [1978]  for  the  situation  where  the  system  of 
partial  differential  equations  had  an  invariant  region  and  by  Hale  [1986]  for 
the  general  situation.  Due  to  the  generality  of  the  methods  in  the  latter  work, 
the  theory  applies  equally  as  well  to  functional  differential  equations  or  delay 
equations  with  diffusion.  The  basic  result  is  that  no  patterns  exist  if  the 

diffusion  coefficients  are  sufficiently  large. 

Once  the  results  arc  known  to  be  valid  for  large  diffusion  coefficients, 

the  next  step  is  to  try  to  understand  the  occurence  of  qualitative  changes  in  the 
flow  through  bifurcations  as  the  diffusion  coefficients  become  smaller.  If  all 
diffusion  coefficients  arc  allowed  to  be  completely  arbitrary,  many  complications 
occur  and  it  is  therefore  natural  as  a  first  step  to  allow  some  diffusion 

coefficients  to  become  small  while  others  remain  very  large.  For  a  system  of 

two  equations  with  special  types  of  nonlinearities  that  occur  in  biological  and 
ecological  models,  such  a  theoretical  investigation  has  been  made  for  the 
bifurcation  and  stability  of  equilibrium  solutions  by  Nishiura  [1981]  and 
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Nishiura  and  Fujii  [1985].  With  the  aid  of  numerical  methods,  the  bifurcation 
of  equilibria  for  the  same  system  for  arbitrary  diffusions  has  been  discussed  by 
Fujii,  Mimura  and  Nishiura  [1982]. 

In  the  work  of  Nishiura  [1982]  and  Nishiura  and  Fujii  [1985],  if  d,,d2 
are  the  diffusion  coefficients,  an  important  role  was  played  by  a  limiting 
system  called  a  "shadow"  system  which  is  obtained  by  letting  d2  — •  ®  and 
consists  of  a  single  reaction-diffusion  equation  with  diffusion  coefficient  dj 
coupled  with  an  ordinary  differential  equation.  It  was  shown  that  the  existence 
and  stability  properties  of  equilibrium  for  the  shadow  system  were  reflected  in 
the  original  pair  of  reaction-diffusion  equations  for  the  diffusion  coefficient 
d2  large.  It  is  the  purpose  of  this  paper  to  carry  this  idea  further  by  showing 
that  the  existence  of  a  compact  attractor  for  the  shadow  system  implies  the 
existence  of  a  compact  attractor  for  the  original  system  if  the  diffusion 
coefficient  d2  is  large. 

Let  us  now  be  more  precise  in  the  statement  of  the  results.  Let  11  be  a 
bounded  domain  in  RN  ,  with  3d  smooth  and  consider  the  system  of 
reaction-diffusion  equations 

3u/3t  «  DjAu  +  f(u,v) 

(Li) 

3v/3t  »  D2Av  +  g(u,v)  in  d 
subject  to  Neumann  boundary  conditions 


(1.2)  3u/3n  -  0  , 

3v/3n 

-  0 

on  3d 

where  u  €  Rm  , 

v  €  Rn 

are 

vectors. 

D,  -  diag(dn,  ,dlm) 

D2  =  diag(d2I,  ,d2n) 

where 

each 

djk  >  °  , 

f  :  Rm  x  Rn  — *  Rm  , 
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g  :  Rm  x  Rn  — *  Rn  are  C2-functions. 

Let  X  «  L2(0,Rm)  ,  Y  «  LJ(n,Rn)  .  Using  the  operator  A  =  -A  ,  one 
can  define  the  usual  fractional  power  spaces  X®  ,  Y®  .  For  convenience  in 

notation,  let  us  assume  that  N  <  3  and  choose  K  <  a  <  1  .  The  latter 

choice  of  a  is  made  to  ensure  that  X®  c  W1,2(n,Rm)  n  L^f^R"1)  , 

Y«  C  Wx'2(n,Rn)  n  L*(n,Rn)  .  For  N  -  1  ,  we  can  take  a  «  K  .  One  can 

then  show  (see,  for  example,  Henry  [1980])  that,  for  any  (u0,vQ)  €  X®  x  Y®  , 
a  >  X  ,  there  is  a  unique  solution  (u(t, -,u0,v0),v(t, -,u0,v0) )  €  X®  x  Y®  of  (1.1), 

(1.2)  through  (uQ,  vQ)  at  t  *■  0  which  is  continuous  in  (t,u0,vQ)  . 

A  set  AC  X®  x  Y®  is  invariant  under  (1.1),  (1.2)  if 

(u(t, -,A),v(t, -,A) )  =  A  for  t  >  0  .  The  set  A  is  a  compact  attractor  if  it  is 
compact,  invariant  and  there  is  a  neighborhood  U  of  A  such  that  the  u>- 
limit  set  of  u  is  A  .  The  u-limit  set  u(U)  of  U  is  defined  as 

u(U)  =  n  cl  u  (u(t, -,U),v(t, -,U))  . 

TfO 

By  the  shadow  system  of  (1.1),  (1.2),  we  mean  the  system 
3u/3t  =  D.Au  +  f(u,{) 

(1.3) 

dt/dt  =  Iftl'1  J  g(u(-,x),Odx  in  n 

n 

with  the  boundary  condition 

(1.4)  3u/3n  ■  0  on  30  . 

The  concept  of  a  compact  attractor  in  X®  x  Rn  for  (1.3),  (1.4)  is  defined  in  a 
manner  analogous  to  the  one  for  (1.1),  (1.2). 
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If  d°  >  0  is  a  given  positive  constant,  we  use  the  notation  Dj  *  d2In  . 
In  general,  we  let  N(S,A)  be  the  5-neighborhood  of  a  set  A  of  a  Banach 


space. 

We  will  impose  the  following  hypothesis: 

(H)  Suppose  there  is  a  compact  set  K  C  X®  *  Rn  and  a  constant  S0  >  0 

such  that  (1.3),  (1.4)  has  a  compact  attractor  An  C  K  , 

i 

u(N(60,ADi  ))  =  Ap^  for  every  Dt  >  d°Im  . 

Theorem  1.  If  (H)  is  satisfied,  then  there  is  a  d2  such  that,  if  D2  *  d°In  , 

then  there  is  a  compact  attractor  An  ,n  C  X“  x  Ya  for  (1.1),  (1.2)  and,  for 

ui  2 

€  >  0  ,  there  is  an  n  >  0  such  that  AD  D  C  N(€,AD  x  {0))  if  D2  >  nln, 

where  {0}  in  AD  x  (0)  means  {0}  C  with  Ya  *  Rn  x  Y^. 

The  last  statement  of  the  theorem  asserts  that  the  attractor  An  n 

1  ,U2 

cannot  be  much  larger  than  AD  x  {0}  if  D2  is  large.  On  the  other  hand, 

it  could  be  smaller  if  we  make  no  further  hypotheses  about  the  flow  on  An 

l 

of  (1.3),  (1.4). 

In  Section  2,  we  prove  Theorem  1.  The  proof  follows  some  of  the  ideas 
in  Hale  [1986]  except  that  a  new  argument  must  replace  the  use  of  Lyapunov 
functions  to  obtain  a  priori  bounds  of  the  solutions  of  the  full  system.  These 
functions  seem  to  be  of  little  use  due  to  the  fact  that  one  cannot  obtain 
information  about  a  perturbed  partial  differential  equation  using  the  derivative 
of  the  Lyapunov  function. 

In  Section  3,  we  discuss  the  difficulties  involved  in  obtaining  An  n  as 

Ul  ’U2 

a  graph  over  (u,z)-space.  Some  restricted  conditions  for  the  existence  of  a 

graph  also  are  given. 


„ '  "  "  A  -v 


Under  the  assumption  that  the  shadow  system  has  a  compact  attractor  with 
certain  properties.  Theorem  1  asserts  the  existence  of  a  family  of  compact 


attractor  AD  ^  for  the  full  system  which  are  upper  semicontinuous  at 
d2  =  •  .  It  is  natural  to  ask  the  opposite  questions:  suppose  that  the  full 

system  has  a  compact  attractor  for  each  D2  *  d°In  .  Does  the  shadow  system 

have  a  compact  attractor  A*  and  is  the  set  {AD  D  |D2  *  djI^A*,}  lower 
semicontinuous  at  d2  =  “  ?  Section  4  is  devoted  to  a  discussion  of  conditions 
which  ensure  that  this  is  the  situation. 

In  Section  5,  we  discuss  the  relationship  between  the  shadow  system  and 
PDE’s  with  nonlocal  spatial  effects  and  hereditary  dependence. 

Theorem  1  can  be  considered  as  a  first  attempt  to  understand  the  behavior 
of  the  solutions  of  systems  of  reaction-diffusion  equations.  Further  information 
could  be  obtained  in  the  following  way.  In  Theorem  1,  there  is  a  restriction 

that  Dj  *  d°Im  .  In  any  particular  problem,  one  could  first  try  to  analyze  the 

shadow  system  for  all  Dx  >  0  and  obtain  a  clear  picture  of  the  dynamics  of 
the  attractor.  Taking  the  limit  as  Dx  — ►  0  gives  some  information  about  the 
types  of  singular  solutions  that  can  occur  at  Dj  =  0  .  Ideally,  one  would  then 
hope  to  obtain  an  attractor  AD  for  D2  *  d°In  and  all  Dj  >  0  .  This 

will  involve  a  very  difficult  analysis  of  the  existence  and  stability  of  large 
amplitude  singular  solutions  near  ■  0  .  These  solutions  will  play  an 

important  role  in  understanding  the  global  flow  for  the  original  equations  in 
the  region  in  (DrD2)  —space  where  D2  is  not  very  large.  Such  a  program 
has  been  partially  carried  out  analytically  and  numerically  for  a  system  of  two 
equations  modeling  problems  in  ecology  by  several  Japanese  mathematicians  (see 


A 
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the  survey  article  of  Nishiura,  Fujii  and  Hosono  [I98hji  T  hc\  ha-,  e 
detail  the  bifurcation  curves  in  (D,,D2)-spacc  for  tnc  one  aril  t« 
equilibrium  solutions  and  have  discussed  the  stability  of  these  solutions 
results  explain  well  the  formation  and  coexistence  of  stable  pattern1 
manner  in  which  the  equilibria  are  dynamically  connected  has  n> 
discussed  and  is  certainly  an  important  factor 

Various  generalizations  of  Theorem  1  arc  possible  For  example,  dilferet.t 
boundary  conditions  may  be  allowed.  Consider  the  equation  tlh  with  tic- 
boundary  conditions 


I  hr; 


bee: 


(1.5) 


Djdu/dn  +  Ej(x)u  =  0 
D,3v/3n  +  E,(x)v  =  0  in  3fl 


where  ErE2  are  diagonal  matrices  and  let 

(1.6)  -  -  |n  r1  J  E2(x)  dx  . 

an 

The  appropriate  shadow  system  for  the  system  (1.1),  (1.5)  with  D2  large 
is  the  system 


(1.7) 


3u/3t  =  DjAu  +  f(u,z)  in  n 

dz/dt  «  l2z  +  inp1  J  g(u,z)dx 

n 


with  the  boundary  conditions 


(18) 


Dj3u/3n  +  Ej(x)u  «  0  in  3fi 


Following  essentially  the  same  proof  as  below  for  Theorem  1  and  the 
estimates  on  e^J  from  Fia le  and  Rocha  [1987a, b],  one  obtains  (X“  ,Y 


o  w 


1 1  a 


Hale  i  ;  ogr  ,  t  r.  a '  the  at’  -  e  ri  _  •  r.  .  !  e..,  a 
Jii:  ere  ntial  cguj!i  n-  w  ,th  4 . :  u  .  r. 


t  i  9i 


3u  3t  -  D,Au  ♦  f‘  u(  .v  , 

d\  b\  -  DjAu  ■*•  g'u(.\ti  in  n 


with  the  boundary.  conditions  M?i  or  (15)  where  the  notation  ur\t  designate 
u  t  >'  8.x  »  =  uu-e.x  I  .  x^S.x)  =  v  ( t  +  0 .  x )  .  0  €  [  -  r  .0 )  .  The  functions  ut-vt  arc 
supposed  to  belong  to  the  spaces  C([-r,0].Xa)  ,  C((-r,0],Ya)  .  The  shadow 

ssstem  is 


3u  3t  =  DjAu  +  f ( u t . z t )  in  ft 

(1.10) 

dz  dt  =  SjZ(t)  +  IDp1  J  g(ut,zt)dx 

n 

i,e  ,  a  functional  partial  differential  equation  coupled  with  a  functional  ordinary 
differential  equation.  We  do  not  state  the  precise  result  on  the  existence  of  an 


attractor  for  (1.9)  which  would  be  analogous  to  Theorem  2. 
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2.  PROOF  OF  THEOREM  1 

Let  Z  C  Ya  be  the  linear  space  of  the  constant  functions, 

Y“  =  Z  «  Y^  ,  v  =  z  +  w  where  z  €  Z  ,  w  €  Y^  , 

(2.1)  z  =  1 0 1 -1  f  v(x)dx  ,  f  w(x)dx  -  0  . 

n  n 

We  can  identify  Z  with  Rn  and  therefore  will  consider  z  as  an  clement 
of  Z  as  well  as  a  vector  in  Rn  . 

If  u(t,  ■)  ,  v(t,  )  are  solutions  of  (1.1),  (1.2)  and  v(t,  •)  =  z(t)  +  w(t,  •)  , 
z(t)  €  Z  ,  w(t,  •)  €  Y^  ,  then 

3u/3t  =  DjAu  +  f(u,z+w) 

(2.2)  dz/dt  =  1 0| _1  J  g(u,z  +  w)dx 

n 

6w/9t  =  D2Aw  +  g(u,z+w)  -  Ifi]'1  J  g(u,z  +  w)dx  in  O 

n 

with 

(2.3)  9u/9n  =  0  ,  3w/9n  =  0  in  9(1  . 

We  are  going  to  consider  this  equation  as  a  perturbation  of  the  system 
3u/3t  =  DjAu  +  f(u,z) 

(2.4)  dz/dt  =  | ni'1  J  g(u,z)dx 

n 

3w/3t  =  D2Aw  in  fl 

with  the  boundary  conditions  (2.3).  By  hypothesis  (H),  (2.3),  (2.4)  has  a  compact 
attractor  An  x  (0)  . 

ui 
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Lct  us  rewrite  (2.2)  as 

3u/3t  *  DtAu  +  f(u,z)  +  F(u,z,w) 

(2.5)  dz/dt  -  \a\-1  J  g(u,z) dx  +  Z(u,z,w) 

n 

3w/3t  =  D2Aw  +  g(u,z)  -  | fl | -1  J  g(u,z)dx  +  G(u,z,w)  . 

O 

Then  using  the  fact  that  a  >  X  and  standard  types  of  estimates  (see,  for 
example,  Hale  [1986]  or  Hale  and  Rocha  [1986a]),  there  is  a  constant  kg  such 
that 

||f(u,z,w)||lJ  ,  |Z(u,z,w)|  ,  |G(u,z,w)|lj  <  ks|!w ]|y-L 

(2.6) 

||g(u,z)  -  |n|-x  J  g(u,z)dx ||  j  <  kB 
0  L 

for  (u,z,w)  €  N(6,K  x  {0})  . 

There  are  jz  >  0  ,  kj  >  0  ,  such  that,  for  d2j  >  d  >  0  ,  we  have 

II'^HIy-L  <  kic'dMt|wlYl  ’  t>0 

1  a  1  a 

(2.7) 

ll'^Hlyl  <  kie‘dMtt'a|w|L2  ,  t  >  0  . 

(see,  for  example.  Hale  [1986]). 

If  we  now  use  the  variation  of  constants  formula  on  the  equation  for  w 
in  (2.5),  then,  as  long  as  (u(t),z(t),w(t) )  €  N(8,AD  x  (0))  ,  we  can  use  the 
relations  (2.6),  (2.7)  to  obtain 


HoLl  *  k1e',£dt  ||w(0)||  .  +  kik6  J  (t-s)'ae'*id(t*s)ds 
a  Y  a  0 

+  J  ||w(s)||  ids 

o  Ya 

Choose  0  <  a  <  n  and  let  n(t)  =  eodt  ||w(t)||  .  ,  y(t)  =  sup{n(s),  0  <  s  <  t) 

Y  a 

and 

CO 

L  =  J  s-<xe*(1-°/#>sds  . 

o 

Then 

n(t)  <  kjC'^"0^*  ti(0)  +  kjkge^^CMd)®-1  +  kjkgL^d^yd) 

Now  choose  d  so  that 

(2.8)  0=1-  kjkgUtid)0'-1  >  0  . 

Then 

y(t)  <  e-1k1e’<'t_o)dt»»(0)  +  0"1  kjkgeodt  Lfpid)**"1  . 

Since  n(t)  -  eodt  ||w(t)||  i  ,  we  have 

*a 

(2-9)  ||w(t)||  .  <  0-lk,e-*idt||w(O)|  .  +  e^kjkgUMd)®'1  . 

1  a  Y  a 

Remember  that  (2.9)  is  valid  for  all  t  €  [0,T]  if  (u(t),z(t),w(t) )  € 

N(6,Ad^  x  (0))  for  t  €  [0,T]  . 

If  (u(t),z(t),w(t))  €  N(6,Ad^  x  (0))  for  all  t  >  0  ,  then  inequality 

(2.9)  implies  that 


I 


(2.10) 


||w(t)||  1  <  e^kj  ||w(0)||  I  +  e^kjkgUftd)®'1  ,  t  >  0 


lim  sup  ^  |]w(t)|j  i  <  G^kjkgLfMd)' 


Therefore,  w(t)  can  be  made  as  small  as  desired  by  taking  w(0)  small  and 
d  large  and  the  lim  sup  of  ||w(t)||  j_  can  be  made  smaI1  by  makin8  d 
large. 

Now  let  us  obtain  a  priori  bounds  on  u(t)  ,  z(t)  ;  namely,  the  solutions 
of 


(2.11) 


3u/6t  =  DxAu  +  f(u,z)  +  F(u,z,w(t)) 


dz/dt  =  |  n  |  -1  J  g(u,z)dx  +  Z(u,z,w (t)) 
ft 


Choose  constants  €,n  ,  0  <  2h  <  €  <  5  ,  and  let  t0(n)  be  the  constant 

such  that  the  solution  (u(t),z(t))  of  (1.3)  (1.4)  with  (u(0),z(0))  G  N(6j,Ad^) 

satisfies  (u(t),z(t))  €  N(n,AD^ )  for  t  >  t0(n)  ,  which  is  ensured  to  exist  by 
(H)  .  There  is  a  constant  <T  such  that,  if  ||F(u,z.w)||  ^  ,  |Z(u,z,w)|  <  T  , 

then  any  solution  (u(t),z(t) )  of  (2.11)  with  (u(0),z(0))  €  N(e,AD^)  must  stay 
in  N(6,AX)  for  0  <  t  <  t0(h)  and  satisfy  (u(t0(h)),z(t0(h)))  €  N(2n,ADj). 
Therefore,  the  same  will  be  true  for  t  €  [t0(n),2t0(n)]  ,  etc.,  and  the  solution 
will  remain  in  N(6,An  )  for  t  >  0  .  To  obtain  this  estimate  on  F,  Z, 

ui 

choose  ||w(0) ||  i  so  small  and  d°  so  large  that  the  right-hand  side  of  the 
first  inequality  in  (2.10)  is  less  than  r/k6  . 

Therefore,  we  have  shown  that  there  is  a  neighborhood  U  of  AD^  x  {0} 
such  that  the  solution  (u(t),z(t),w(t))  of  (2.2),  (2.3)  with  initial  data  in  U 


mmm 

mm 

s 

mm 

ss 
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and  dj  sufficiently  large  stays  in  NCS.Aj  x  {0})  for  t  H  .  Thus,  u(U) 
is  compact  (see  Henry  [1980],  Th.  3.3.6)  and  uj(U)  =  AD  is  a  compact 

attractor  for  all  D,  >  d?Im  and  D,  >  d?L  . 

1  1  m  2  2  n 

From  the  second  inequality  in  (2.10),  we  can  use  an  argument  similar  to  the 
one  above  to  show  that,  for  any  sequence  — *  0  as  j  — *  •  ,  there  is  a 

sequence  d2-  — *  ®  as  j  — *  •  such  that 

A  =  U(U)  €  N(n  Ad  x  {0})  . 

d1d2  j  ui 

This  implies  the  last  assertion  in  Theorem  1  and  the  proof  is  complete. 


3.  THE  ATTRACTOR  AS  A  GRAPH 


In  this  section,  we  discuss  the  possibility  of  the  attractor  AD  D 

being  a  graph  over  (u.z)-space.  It  is  tempting  to  try  to  obtain  such  a  graph  by 

applying  the  method  of  center  manifolds  to  the  system  (2.2),  (2.3)  attempting  to 

obtain  an  invariant  manifold  which  contains  An  n  and  is  defined  bv 

U1’U2 

w  =  h(u,z)  for  some  function  h  .  After  appropriately  re-defining  f  ,  g 

outside  a  neighborhood  of  the  attractor  so  that  they  are  bounded  functions,  such 
an  integral  manifold  would  be  required  to  satisfy  the  equations 


9u/9t  *  DjAu  +  f (u,z  +  h(u,{)) 


dz/dt  *  |n|'1  J  g(u,z  +  h(u,U)dx 
n 

u(0)  =  u0  ,  z(0)  =  z0 


aei  0 

h(u0,z0)  =  (Th)(u0,z0)  s  J 


d«f  -D.As 


g(u(s),z(s)  +  h(u(s),z(s)) )  - 


W1  J  g(u(s,x),z(s),h(u{s,x),z(s)))dx  jds 


To  define  this  integral  operator  Th  requires  that  u(s)  must  be  defined  on 
(-°°,0]  .  One  does  not  expect  such  solutions  to  exist  for  all  values  of  uQ 
because  of  the  smoothing  properties  of  the  solutions  of  9u/9t  =  DjAu  . 

These  remarks  seem  to  indicate  that  the  standard  method  of  center 
manifolds  will  not  apply  to  this  problem. 

Another  approach  that  could  make  use  of  center  manifold  techniques  is  to 
assume  that  An  for  the  extended  f  ,  g  lies  in  a  finite  dimensional 

ui 

invariant  subspace  which  is  normally  hyperbolic  for  the  shadow  system.  This 


WWW 
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can  be  accomplished,  for  example,  if  the  spectrum  of  the  Laplacian  has 
sufficiently  large  gaps.  The  proof  can  be  supplied  as  in  Mallet-Paret  and  Sell 
[1986]  for  inertial  manifolds  or  one  can  use  the  integral  equation  approach  for 
center  manifolds  in  a  form  similar  to  that  mentioned  above.  This  implies,  for 
example,  that  the  attractor  will  be  a  graph  if  the  equation  for  u  acts  in 
only  one  space  variable. 


4.  LOWER  SEMICONTINUITY  OF  THE  ATTRACTORS 


In  this  section,  we  assume  the  existence  of  compact  attractors  An  n  for 

u  j 

(1.1),  (1.2)  satisfying  certain  properties  and  conclude  that  the  shadow  system 
(1.3),  (1.4)  has  a  compact  attractor. 


Lemma  1.  For  t  >  0  .  T  >  0  ,  there  exists  d2  *  d2(£,T)  such  that:  if 

(u(t),v(t))  and  (w (t),$(t))  are  solutions  of  (1.1).  (1.2)  and  (1.3).  (1.4)  respectively 
with  u(0)  =  w(0)  .  v(0)  =  ((0)  then,  for  D2  >  d2(e,T)In  ,  the  following  are 

valid: 

sup  I  u(  t )  -  w(t)|  <  e  ,  sup  |v(t)  -  ((t)|  <  £ 

t€(0,T]  x  t^O.T) 

sup  I  V  ( t )  -  V(t)|  <  £ 

t^O.T)1  'YI 

where  V(t)  ■  1 0 1 " 1  J  v(t)dx  . 

n 


Proof.  By  gj ))  ,  4>  €  X®  ,  0  €  Ya  ,  we  denote  the  projection  of  the 

function  g(0(  ),i l>{  ))  onto  Yj  along  Z  .  Let  us  first  give  an  estimate  on 
vjjt)  =  v(t)  -  V(t)  ,  which  satisfied  the  following 


D2A(t-s) 

e 


g^(u(s),v(s)  +  v_L(s))ds 
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lv±(t>||Yj  «  *i  /(t-s)-«c-^^)||gl(u(s),v(s))||YOds 

(4.1) 

t  kjMm-aVC^dz)1*®  for  t  €  [0,T] 

where  M  »  sup{|Jgx( u(s),v(s) )  (Jy0  j  0  <  s  <  r}  . 

Now  rewrite  the  equations  for  ( u,V )  as  follows: 

3u/9t  =  D,Au  +  f(u,V)  +  Rj(v^) 

dV/dt  -  inp1  j  g(u(y,t),V(t»dy  +  R2(vj^) 
n 

where  Ri(vl>  -  f(u,V+V-L)  -  f(u,v) 

and 

R2(vl)  ■  I0!'1/  [«(u(y,t),V(t)+vj_(y,t))  -  g(u(y,t),v(t))]  dy  . 

n 

By  using  the  estimate  in  (4.1),  one  can  find  a  constant  K  >  0  such  that 
sup  ||Ri(v±)|[  <  Kdf’1  ,  i  -  1,2  . 

t^0,T]  "Y 

These  estimates  prove  the  existence  of  d,  ■  d2(c,T)  in  the  lemma. 

Let  us  denote  by  the  semiflow  on  Xa  x  Ya  generated  by 

(LI),  (1.2)  and  make  the  following  hypotheses: 

M*l)  Td  d  (t)  has  a  compact  attractor  An  n  *  <P  in  X®  x  Y“  , 

1’  2  ul'u2 

(H-2)  For  any  bounded  set  B  C  X“  x  Ya  and  any  €  >  0  ,  there  is  a 
lo  *  t0(c,B)  such  that: 

tDi.dO)B  C  n(«,aDiD!)  ,  V,  >,0 

(H-3)  Cl  U  Ad  d  'S  compact  for  some  d2  >  0  . 

Dj»d°In  v  2 
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Theorem  3.  Under  the  hypotheses  (H-l)-(H-3),  the  shadow  system  (1.3),  (1.4)  has  a 

compact  attractor  An  and ,  moreover, 

*^1  00 
*  I 

Ajj  =  O  Cl  U  Ap  p 

6)dj  Dj>5ln  Dl'°2 


Proof.  n  Cl  u  An  n  is  non-empty  and  compact  because  of  (H-l)and  (H-3). 

8)dj  Dj«in  11  2 

Let  us  denote  by  Tn  (t)  the  semiflow  on  X®  x  Z  generated  by  the  shadow 

1 1 

system  (1.3),  (1.4). 

Firstly,  we  show  that  Tn  (t)  has  a  compact  attractor  An  and  that 

U\  CO  Ul  CO 

a 

Ad  C  o  Cl  U  Ad  D  holds.  For  any  bounded  set  B  C  X  x  Z  ,  the 
6>d°  d,*8i  r  2 

2  2  n 

Set  {TD1,D2(t)B|t  >  0}  is  bounded  by  virtue  of  (H-l).  We  therefore  repeat  the 
arguments  in  the  proof  of  Lemma  1  uniformly  in  the  initial  data  (<(>,$)  €  B  , 
namely,  for  any  e  >  0  and  T  >  0  ,  there  exists  a  constant  60  =  80(£,t,B) 
such  that 


dis,xaxYa  ^TD  j  <  € 

for  t  €  [0,T]  ,  D2  *  60(T,e,B)In  .  By  the  hypothesis  (H-2),  there  is  a  t0(e,B) 
such  that  TDi  D^(t)B  C  n[«,Ad^d  ]  for  t  *  t0(«,B)  D2  >  d°In  .  For  any 
increasing  sequence  {T.)'l1 ,  T.T®  ,  jt®,  let  j0  =  j0(£)  be  the  least  index 
j  for  which  T.  >  t0(€,B)  and 

n[.,  A  Cl  U  A  1  D  Cl  U  {a  ID,  >  60«,T  B)l, ) 

I.  6>d?  D,>81  1  lJ  A  13  ' 

4  2  n 


s  v o  ■o'  o*o  ■C**1 va'aV  ^ ^ 
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For  this  choice  of  j0(€)  ,  we  have: 

(1)  dist(TDiDj(Tj^€j)B,TD^  J  Tjo(£))B)  <  e  for  D2  >  6o(c’’rj0(€)'B)In 

(2)  td1,dj(tj0(€))b  c  n(€.ad1iDj)  .  VD2  >  <*?.  . 

(3)  N(c,  n  Cl  U  AD  D  )  D  Cl  U  {AD  |D2  >  80(«,t.(6),B)Ib}  • 

6>dj  Dj)8ln  12  12  0 

From  these  three  properties,  it  follows  that 


rD  (r  fe))B  c  N(3«,  n  cl  u  ad  D  ) 
J°( °  8>d«  d,«i  r  2 

2  2  n 


Let  jk  =  j0(l/k)  ,  k  >  1  .  Hence,  the  limit 


lim 

k-*«o 


11  Td  JTj^B  C  n  0CI  U  Adi'D 

-  *•  k  6>d°  D2>8ln  1 


exists. 


Since  the  sequence  (Tj}j°=i  >  TjT®  >  jt«  is  arbitrary,  we  have  thus 
established 


(4.2) 


UT  (B)  -  r,ClUTD  (s)BC  n  Cl  U  A 
Dl,0°  t*0  S*t  8?d,  D2^5ln 


This,  in  particular,  implies  that  Tn  *(t)  is  bounded  dissipative  and 

ui • 

{Td  *(t)B|t  >  0}  is  bounded  for  B  C  Xa  *  Z  bounded.  One  C3n  also  show 
that  Tn  oo( t )  is  compact  for  t  >  0  (see  Hale  [1985]  for  detail).  These  three 

l  ■ 

properties  are  sufficient  to  guarantee  that  T^  „(t)  has  a  compact  attractor 
An  .  (see  Hale  [1985]).  From  (4.2),  it  follows  that 
AD  .  c  n  Cl  u  AD  D  . 

6>d j  Dj>6Jn  r  2 


In  order  to  prove  A 


,  it  suffices  to  note  that 


Di  "  3  °  oC*  U  A°i'D2 

1  6>d°  d,>6i  1  2 

2  2  n 


Lemma  1  actually  shows  that  the  semiflow  Tn  n  (t) 

ul'u2 

infinity.  This  fact  together  with  the  existence  of  An 

ui 

to  ensure  that  An  n  is  uppcr-semicontinuous  in  D, 

ul’u2  1 


A D  .  I>  n  Cl  U  ad  d  . 

11  6>d°  D2>6ln  r  2 


This  completes  the  proof  of  Theorem  3. 


is  continuous 
»  and  (H-3) 
at  infinity. 


in  D2  at 
is  sufficient 
namely, 
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5.  SHADOW  SYSTEMS  AND  FUNCTIONAL  DEPENDENCE 


For  a  system  of  two  reaction-diffusion  equations  with  one  diffusion 
coefficient  large  relative  to  the  other,  we  have  seen  that  the  flow  on  the 
attractor  can  be  reduced  to  the  discussion  of  the  shadow  system  consisting  of  a 
PDE  coupled  with  an  ODE  with  nonlocal  terms.  In  some  situations,  the  flow  of 
the  shadow'  system  can  be  reduced  to  the  discussion  of  a  scalar  PDE  with 
nonlocal  effects  in  the  spatial  variables.  In  this  section,  we  give  some 
illustrations  of  this  fact. 

Suppose  the  shadow  system  is  given  by 


u  t  =  dAu  +  f(u,0 

=  l^r1  J  g(u(  )  dx  in  (1 
n 


with  the  boundary  condition 


3u/3n  =  0  in  3fi 


Suppose  that  (5.1),  (5.2)  has  a  compact  attractor  Ad  and  that 


g(u,v)  =  -X[v-h(u)] 


where  X  >  0  is  a  positive  constant  and  h(u)  is  a  C2-function. 

If  (u(t)T(t))  €  Aj  for  t  e  R  ,  then,  in  particular,  {(t)  is  a  solution 
of  the  equation 

-  -H(t)  +  X  1 0 ( _1  J  h(u(t,x))dx 
n 

which  is  bounded  on  R  .  Therefore,  $(t)  is  uniquely  defined  by  the  formula 


(5.4) 


H[u](s)ds 


l(t)  =  x  j  e'X(t's) 


where 

(5.5)  H[u](t)  =  inf1  J  h(u(t,x))dx 

n 


This  implies  that  u(t) 
(5.6)  ut  =  dAu  + 


is  a  solution  of  the  equation 
f  0 

f  u(t,),vj  e  H[u](t  +  s)  ds 

v  _  OO 


in  n 


with  the  boundary  condition  (5.2).  This  is  a  retarded  PDE  with  nonlocal  effects 
in  the  spatial  terms. 

In  summary,  if  g  satisfied  (5.3),  then  the  flow  on  the  attractor  for  the 
shadow  system  (5  1),  (5.2)  is  equivalent  to  discussing  the  flow  on  the  attractor 
for  (5.6),  (5.2)  making  use  of  (5.4).  Of  course,  we  must  make  certain  that  (5.6), 
(5.2)  defines  a  semigroup  in  some  appropriate  space. 

To  obtain  a  solution  of  (5.6),  one  must  specify  a  function  9:(-°°,0]  — *  X“ 
and  then  attempt  to  use  (5.6)  to  extend  <f  to  a  function  u(t,<p)  defined  on 
(.<*>,«>)  with  u(t,<p)  =  <p(t)  for  t  <  0  .  There  are  several  natural  spaces  for 
the  <*>  ’s;  for  example,  for  any  0  <  y  <  V  ,  the  space 


C?  =  (<f  €  C ((-®,0),Xa) ;<?(6) e?®  — *  limit  as  9  — »  -°°) 


with  the  sup  norm 


or  the  space 


L7p((-«0),Xa}  x  Xa 

Ly  p((-"D,0),Xa)  *  |<f:(-°>,0)  — •  Xa?9  measurable  and 
0  3 

J  e^p6 1  if(0)  |  pd0  <  « 

.oo  J 
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with  the  norm 


|(W|7,p  =  [|“o!P  +  !°fpQ  |‘P(0)|pde]] 


With  either  of  these  spaces,  one  can  follow  the  usual  procedures  to  obtain 
the  local  existence  and  uniqueness  of  solutions  of  the  initial  value  problem  for 
(5.6),  (5.2). 

Let  u-;  restrict  our  discussion  to  C?  .  If  Ad  C  Xa  x  R  is  a  compact 
attractor  for  (5.1),  (5.2),  then  (5.6),  (5.2)  has  a  compact  attractor  Ad  c  C?  and 
(u(t),((t))  €  Ad  if  and  only  if  \(l)  is  given  by  (5.4)  and  T(t)u  €  Ad  where 
r(t)u(0)  =  u(t+0)  ,  -00  <  0  <  0  . 

It  is  clear  that  properties  of  the  solutions  of  equations  of  the  form  (5.6) 
need  to  be  investigated  in  more  detail.  In  this  paper,  we  are  content  with  a 
few  remarks. 

If  we  let  ft  =  (0,7t)  and  let  uQ  be  the  equilibrium  point  of  (5.6),  then 
the  eigenvalues  of  the  linear  variational  equation  about  u0  are 


Xn  =  a  -  n2d  ,  n  >  1 


XQ  =  a  +  e(l  +  X0)‘1 


with  corresponding  eigenfunctions  e  n  cos  nx  ,  0  €  (-“,0]  ,  x  €  [0,77]  where 

«  =  fu(uo’*o)  »  8  “  fv(uo’*o)  »  *0  =  h(u0)  •  Now  suppose  that 


0  <  o  <  1  ,  B  <  0  ,  a  +  B  <  0  . 


These  inequalities  are  the  usual  ones  corresponding  to  the  Turing 
conditions  for  the  destablization  of  equilibria  in  the  original  pair  of  reaction 
diffusion  equations  (see,  for  example,  Nishiura  [1982]). 


vmmm 


’  >  ,V  •  u'V  V-  v ' '  .■  1 
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If  (5.7)  is  satisfied,  then  ReXQ  <  0  .  If  we  consider  d  as  a  parameter, 
then  each  Xn  <  0  if  n  is  large.  If  we  decrease  d  ,  then  there  is  a 

bifurcation  at  d  =  a  and  another  equilibrium  will  arise  which  is  spatially 

nonhomogeneous.  This  seems  to  be  typical  of  the  Turing  mechanism.  The 

nonlocal  spatial  effects  make  the  eigenfunction  corresponding  to  the  dominant 

eigenvalue  spatially  dependent. 

Another  interesting  equation  is  obtained  by  taking  a  limiting  situation  in 
(5.4).  If  X  — »  *  ,  then  the  equation  (5.6)  should  have  the  dynamics  given  by 
the  simpler  equation 

(5.9)  ut  =  dAu  +  f|u,|fl|_1  J  h(u(t,x))dx 

n 

with  the  boundary  condition  (5.2).  Equations  of  this  type  have  been 

encountered  by  Chafee  [1981],  Levin  and  Segal  [1982],  [1985],  where  they  also 

observed  that  stable  patterns  could  be  generated  by  nonlocal  spatial  effects 
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